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ABSTRACT 

The NASA Dryden Flight Research Center and the 
Technical University of Munich are cooperating in a 
research program to assess the impact of unstable long- 
term dynamics on the flying qualities of planes in hyper- 
sonic flight. These flying qualities issues are being investi- 
gated with a dedicated flight simulator for hypersonic 
vehicles located at NASA Dryden. Several NASA 
research pilots have flown the simulator through well- 
defined steady-level turns with varying phugoid and height 
mode instabilities. The data collected include pilot ratings 
and comments, performance measurements, and pilot 
workload measurements. The results presented in this 
paper include design guidelines for height and phugoid 
mode instabilities, an evaluation of the tapping method 
used to measure pilot workload, a discussion of techniques 
developed by the pilots to control large instabilities, and a 
discussion of how flying qualities of unstable long-term 
dynamics influence control power design requirements. 

INTRODUCTION 

Much of the work in recent years focused on proposals 
and evaluations of hypersonic vehicle concepts. 1 Many of 
these concepts required development of technology in 
fields such as structures and propulsion, and researchers 
have appropriately focused in these areas. However, 
deficiencies also exist in the flying qualities criteria for 
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hypersonic vehicles. One example where the criteria are 
deficient is in the impact of long-term dynamics on hyper- 
sonic flying qualities. 2 The long-term dynamics consist of 
the phugoid and height modes of motion. 

The current military standard 3 does not adequately 
address the phugoid and height mode flying qualities 
requirements at hypersonic speeds. The phugoid flying 
qualities criteria are based on data from the subsonic 
regime. However, substantial differences exist between 
phugoid dynamics in subsonic and hypersonic flight. From 
the constant energy equation, a phugoid mode excited at 
high speeds produces greater altitude excursions for a 
given velocity excursion than a phugoid mode excited at 
low speeds. To make matters worse, tight control of alti- 
tude becomes more important in hypersonic flight than in 
subsonic flight because of engine performance sensitivity 
to dynamic pressure. Thus, the current criteria and data 
collected on the phugoid mode may not be applicable for 
the supersonic and hypersonic speed regimes. 

The height mode is a mode of motion that becomes 
more significant in hypersonic (and supersonic) flight than 
in subsonic flight. The height mode, not addressed by the 
military standard, is typically a first-order response result- 
ing from thrust variation with density gradient or speed. A 
more detailed description of the phugoid and height mode 
characteristics has previously been given. 4 

Berry addressed the issue of long-term dynamics 5 by 
examining the impact of the propulsion system on the 
phugoid and height mode dynamics of a hypersonic lifting 
vehicle developed by Etkin. 6 Berry discovered that air- 
breathing propulsion systems can produce flight regions of 
unstable height and phugoid modes. Berry thought that 
safe pilot control of these instabilities would have strong 
implications on redundancy management and backup sys- 
tems required for the flight control system of hypersonic 
vehicles. The flying qualities impact of these unstable 
modes was, however, beyond the scope of his work. 



The need for this data led the NASA Dryden Flight 
Research Center and the Technical University of Munich 
into a cooperative research program to assess the flying 
qualities impact of unstable long-term dynamics on hyper- 
sonic planes. The primary objective of the research is to 
develop and validate flying qualities criteria for these 
dynamic characteristics of hypersonic vehicles. These 
issues are investigated with a dedicated flight simulator 
for hypersonic vehicles and an SR-71 aircraft located at 
NASA Dryden. Previous results from this cooperative 
research have been documented. 7 "’ 10 The results presented 
in this paper are as follows: 

• the design guidelines for height and phugoid mode 
instabilities 

• the piloting techniques developed to control high 
degree of instabilities 

• the evaluation of the tapping method used to 
measure pilot workload 

• the influence of flying qualities of unstable long- 
term dynamics on control power design 
requirements 

NOMENCLATURE 

A state matrix 

B control matrix 

BTL basic tapping level 

C output matrix 

C D change in the nondimensional drag coeffi- 

cient with respect to speed changes 

C m change in the nondimensional moment coef- 

ficient with respect to flightpath change 

C change in the nondimensional moment coef- 

« ficient with respect to elevator change 

GHAME Generic Hypersonic Aerodynamic Model 
Example 

H altitude, ft 

j a given run at each instability level 

KEAS equivalent airspeed, knots 

K g control system gain, C m /C m ^ 

LTL loaded tapping level 

N total number of runs flown at each instability 


NASP National Aerospace Plane 

PML perceptual motor load 

p roll rate, deg/sec 

P' rms variation of the parameter altitude or 

rms r 

KEAS for a given run 

P rms variation of the parameter altitude or 

rms r 

KEAS across the total number of runs at a 
given instability 

q pitch rate, deg/sec 

r yaw rate, deg/sec 

rms root mean square 

s Laplace operator 

u control input vector 

V true airspeed, ft/sec 

X element in matrix A modified for height 

mode variation 

x state vector 

x change of x with respect to time 

y output vector 

a angle of attack, deg 

p sideslip angle, deg 

y flightpath angle, deg 

A change in a parameter 

& a aileron position, rad 

8 lateral stick deflection, in. 

ap 

elevator position, rad 

8 longitudinal stick deflection, in. 

5 r rudder position, rad 

8 rudder pedal position, in. 

Sj. throttle position, percent of full throw 

C >dr dutch roll damping 

0 pitch attitude, deg 

$ bank angle, deg 
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\y heading angle, deg 

c o dr dutch roll frequency, rad/sec 

SIMULATOR DESCRIPTION 

The space shuttle and the National Aerospace Plane 
(NASP) programs 11 successfully applied the fixed-base 
flight simulator used in this research in previous hyper- 
sonic flight research. For this study, the hypersonic vehicle 
aerodynamics and propulsion system were linearized 


models based on the Generic Hypersonic Aerodynamic 
Model Example (GHAME) at Mach 10 and an altitude of 
110,000 ft. The GHAME model represents a generic, 
unclassified hypersonic vehicle incorporating six degrees 
of freedom; oblate earth equations of motion; and a turbo- 
jet, ramjet, and scramjet propulsion system combination. 
A simple, conventional control system (pitch rate feed- 
back) designed for the linearized model provided Level 1 
longitudinal short-term and lateral-directional flying qual- 
ities at the flight condition (fig. I). 3 The pilot- vehicle 
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(a) Longitudinal short-period flying qualities design 
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(b) Lateral-directional dutch roll flying qualities design. 
Figure 1. Baseline Mach 10 flying qualities designs. 
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(c) Lateral-<lirectional roll mode and dutch roll mode flying qualities design. 

Figure 1. Concluded. 


interface was a conventional center stick and rudder pedal. 
The appendix details the control system and linear equa- 
tions of motion for the flight condition used in this study. 

The cockpit instrumentation was based on a configura- 
tion used for the space shuttle. A central control panel and 
a simulated head-up display (fig. 2) provide the pilot with 
cues adequate to perform hypersonic maneuvering. 
Included in the head-up display were vertical and longitu- 
dinal accelerations, which proved useful for flightpath and 
speed control during space shuttle and NASP research. A 
fixed scale and arrow with rolling digits displayed vertical 
speed while a moving tape and arrow with rolling digits 
inside indicated altitude. A moving scale and fixed arrow 
displayed the heading with a resolution of 0.5°. Because 
of experience gained from initial simulator runs, it was 
decided to display an inertial vertical speed indicator in 
the form of a diamond on the pitch ladder. 

TEST PROCEDURE 

Pilots flew a well-defined steady-level turn to investi- 
gate the flying qualities implications and characteristics of 
various levels-of-height and phugoid instabilities. The 
definition of the maneuver was to perform a 30° bank turn 
and capture a 12° heading change. This maneuver had to 
be performed while immediately eliminating a vertical 
speed initial condition and maintaining constant airspeed 
and altitude. 

Adequate performance was defined as capturing the 
final heading within ± 1°, maintaining speed within 


± 10 knots of the initial speed, and maintaining altitude 
within ± 600 ft of the initial altitude. 

Desired performance was defined as capturing the final 
heading within ± 0.5°, maintaining speed within ± 5 knots 
of the initial speed, and maintaining altitude within ± 300 
ft of the initial altitude. 

These performance criteria were chosen arbitrarily as 
reasonable but tight constraints and are shown in Table 1 . 
The turning portion of the maneuver lasted approximately 
2 min and termination of the maneuver occurred 1 min 
after capturing the final heading. Initial conditions were 
wings level at Mach 10 and an altitude of 110,000 ft. A 
vertical speed initial condition of 20 ft/sec was introduced 
to increase pilot workload. 


Table 1. Adequate and desired performances. 


Steady Level Turn 

Controlled 

Adequate 

Desired 

parameters 

performance 

performance 

Target heading 

±1° 

±0.5° 

Target altitude 

±600 ft 

±300 ft 

Trim speed ( KEAS) 

±10 kn 

±5 kn 


The pilots involved in the test program were experi- 
enced in high-speed flying. Their experience includes a 
wide spectrum of supersonic and hypersonic vehicles, 
including the SR-71, YF-12, and X-15. 
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Figure 2. Simulated head-up display with modifications for hypersonic flight. 


The pilots evaluated phugoid mode instabilities with decides if simple side tasks (for example, radio communi- 

time-to-double-amplitudes down to 1.7 sec and height cations) or complex side tasks (for example, flying in a 

mode instabilities down to 3.5 sec, which were empiri- complex air traffic control environment) are possible, 

cally determined limits of controllability for this study. After these distinctions, the pilot must estimate the level of 

Similar levels of instability were previously studied 12 compensation required to perform the maneuver, resulting 

using a large supersonic transport model with an aperiodic in a rating from 1 to 13. 

phugoid time-to-double-amplitude of 4.6 sec in the 

approach and landing condition. After each run, the pilot A ^ me history recorded from each test run provided 
gave Cooper-Harper ratings 13 (fig. 3), comments, and information on pilot performance. Pilot performance was 
estimations of workload. measured throughout the entire maneuver as the root- 

mean-square (rms) variation from the initial condition of 
Figure 4 shows the workload scale developed for this the key parameters altitude and airspeed. Pilot perfor- 

study to provide insight on pilot workload independent of mance was measured for only those test runs where the 

performance. For this workload scale, the pilot estimates pilot maintained control. The following equation describes 

the complexity of the side tasks possible while performing the summary of pilot performance for all the pilots at each 

the maneuver. If the airplane is controllable, the pilot instability: 
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Aircraft Demands on the pilot in selected 

characteristics task or required operation 



Excellent - Highly 
desirable 

Pilot compensation not a factor for 
desired performance 

„ Good • Negligible 
deficiencies 

Pilot compensation not a factor for 
desired performance 

Fair - Some mildly 

unpleasant 

deficiencies 

Minimal pilot compensation 
required for desired performance 


Is It satisfactory 
without 
improvement? 


Deficiencies 

warrant 

Improvement 


Minor but annoying 
deficiencies 

Moderately objection- 
able deficiencies 

Very objectionable 
but tolerable 
deficiencies 


Desired performance requires 
moderate pilot compensation 

Adequate performance requires 
considerable pilot compensation 

Adequate performance requires 
extensive pilot compensation 


r Is adequate \ 
performance attainable' 
with a tolerable pilot j 
V workload? S 


Deficiencies 

require 

Improvement 


Major deficiencies 

Major deficiencies 
Major deficiencies 


Adequate performance not obtainable 
with maximum tolerable pilot compen- 
sation; controllability not In question 

Considerable pilot compensation 

Is required for control 

Intense pilot compensation is required 
to retain control 


Is It controllable? 


Improvement 

mandatory 


Major deficiencies 


Control will be tost during some 
portion of required operation 


Pilot decisions 


Figure 3. Cooper-Harper rating scale. 
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Complex side tasks (e.g. 


Compensation not noticeable 


flying In a complex ATC 


Minimal but noticeable compensation 


environment without 
any support) and 
flying requires: 


Small but more noticeable compensation 



Moderate compensation 



Considerable compensation 





Simple side tasks (e.g. 


Moderate compensation 


radio communications) 


Considerable compensation 


and flying requires: 


Immense compensation 





No more side tasks 


Moderate compensation 

□ 

possible (backseater 


Considerable compensation 

0 

required) 


Immense compensation 

0 



Not tolerable compensation 

0 




Aircraft runs out of control 

0 


Pilot decisions 


Figure 4. Workload rating scale. 


N 



P' 2 

rms 

L=_i 


1 N 


where 

P = altitude or airspeed 

P' ms = 11115 variation of the parameter about its 

initial condition for each run 
N = the total number of runs flown at each 
instability level 

Pfms = Tms variation about the initial condition 
across all runs at each instability 

The P ^ parameter provides a summary of the perfor- 
mance achieved across all the runs at a given instability. 


The phugoid and height mode instabilities were imple- 
mented through unstable feedback in the control system 
and modification of the aerodynamic models. Phugoid 
mode instabilities were implemented by positive feedback 
of flightpath angle (fig. 5). Height mode instabilities were 
implemented by modifying the drag coefficient due to 
speed, C D (fig. 6). This implementation allowed phugoid 
and height mode instabilities to be set independently from 
each other. The eigenvectors of the phugoid and height 
modes implemented with these methods were verified to 
have accurate characteristics. All phugoid instabilities 
evaluated were in the aperiodic region; however, the 
baseline phugoid case was periodic and neutrally stable. 
The baseline height mode case was stable with a 
time-to-half- amplitude of 77 sec. Tables 2 and 3 show the 
phugoid and height mode instabilities evaluated. 
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Table 2. Phugoid mode instabilities evaluated. 



Root location, 
rad/sec 

Time to double 
amplitude, sec 

0 

0 

— 

4.2 

0.05 

13.6 

10 

0.16 

4.4 

20 

0.30 

2.3 

25 

0.36 

1.9 

30 

0.42 

1.7 


Table 3. Height mode instabilities evaluated. 


V 

Root location, 
rad/sec 

Time to double 
amplitude, sec 

0.00001 

-om 

— 

-0.00003 

0.02 

40.8 

-0.00006 

0.04 

18.2 

-0.00008 

0.05 

13.6 

-0.00018 

0.12 

5.9 

-0.00024 

0.16 

4.4 

-0.00030 

0.20 

3.5 



Figure 5. Aperiodic phugoid mode instability characteristics. 
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Figure 6. Aperiodic height mode instability root contours 







To increase pilot workload, the test runs included atmo- 
spheric density perturbations for the pilot to fly through. 
The disturbances occurred within 1 to 20 sec of each other 
in the shape of a cosine wave lasting 20 sec. The amplitude 
of the disturbance, whose sign was randomly modelled, was 
2.5 percent of the standard density at an altitude of 110,000 
ft and consistent with atmospheric variations observed from 
YF- 12 flights. 14 

During some of the test runs, a secondary side task for the 
pilot was instituted to measure pilot workload. The method, 
called tapping, required the pilot to tap a button located on 
the stick as rhythmically as possible. The theory states that 
as the pilot workload increases, the irregularity of the tap- 
ping increases. A comparison between the regularity of the 
tapping for a given run (loaded tapping level, LTL) and the 
regularity of the tapping during straight and level flight with 
no disturbances (basic tapping level, BTL) yields a measure 
of the pilot workload or perceptual motor load (PML). The 
PML is calculated as follows: 

PML = (LTL -BTL)/BTL 

Information on the tapping method as it was applied to 
this research has previously been given. 10 General informa- 
tion on the tapping theory has also been reported 
previously. 15 

RESULTS AND DISCUSSION 

The following section discusses how phugoid and height 
mode instabilities affect Cooper-Harper ratings, compares 
the effects of workload and pilot performance on Cooper- 
Harper ratings, and analyzes pilot technique in controlling 
large instabilities. The effect of control power on 


Cooper-Harper ratings is examined, and the tapping 
method is evaluated. 

Phugoid and Height Instability Impact on Cooper- 
Harper Ratings 

Figure 7 shows the Cooper-Harper ratings as a function 
of phugoid real-root locations without use of the tapping 
method. Five data points, generated by three pilots, were 
averaged at each instability. Eighty-three percent of the 
ratings fell within a ± 1 pilot-rating band, illustrating the 
consistency of the results. The ratings stay relatively con- 
stant at low phugoid instabilities until a root location of 
0.15 rad/sec, where the ratings degrade sharply. Level 2/3 
and Level 3/Uncontrollable borders occur at unstable root 
locations of 0.23 and 0.44 rad/sec respectively, or 3 and 
1.4 sec time-to-double-amplitude respectively. 

Figure 8 shows the nontapping Cooper-Harper ratings 
as a function of height mode instability. Seven data 
points, generated by four pilots, were averaged at each 
instability. Seventy-three percent of the ratings fell within 
a ± 1 pilot-rating band, illustrating the consistency of the 
results. Similar to the phugoid data (fig. 7), relatively 
constant ratings occur for the low instabilities, degrading 
sharply as instability increases above a root location of 
approximately 0.015 rad/sec. This result indicates a 
“cliff-like” characteristic in the height mode data. Com- 
paring the slopes of the height mode and phugoid mode 
Cooper-Harper data (figs. 7 and 8) indicates that the 
“cliff” characteristic exists in the phugoid data to a lesser 
degree. Level 2/3 and Level 3/Uncontrollable borders 
occur at unstable root locations of approximately 0.045 
and 0.14 rad/sec respectively, corresponding to 15 and 5 
sec time-to-double-amplitude. 


Cooper-Harper 

rating 



Figure 7. Average Cooper-Harper rating as a function of phugoid instability (nontapping data). 
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Figure 8. Average Cooper-Harper rating as a function of height mode instability (nontapping data). 


Assuming Level 2 ratings are acceptable for backup 
control modes of hypersonic vehicles, these data indicate a 
time-to-double-amplitude of 15 sec for the height mode 
and 3 sec for the phugoid mode are acceptable limits. 
However, the apparent “cliff-like’ 5 tendencies in the rat- 
ings may warrant extra margin in design. Comparison of 
the phugoid and height mode data suggests the pilot is 
more sensitive to the height mode instability than the 
phugoid mode instability, as the height mode borders 
occur at larger time-to-double-amplitudes. This result 
might be expected because the height mode produces a 
divergence in energy as well as flightpath, thus requiring 
both throttle and pitch control to correct a deviation. The 
phugoid mode, on the other hand, is primarily a constant 
energy mode and deviation can be corrected with the pitch 
controller alone. 

Comparison of Workload and Performance Effects on 
Cooper-Harper Ratings 

In the previous section, all of the evaluations included 
both controllable and uncontrollable cases. In this section, 
the comparison of performance measures and correspond- 
ing pilot ratings includes only the controllable cases 
because the performance was calculated only for those 
cases. 

The trends in the Cooper-Harper rating data and the 
workload rating data for the controllable phugoid instabil- 
ities without tapping (fig. 9) are similar. Both Cooper- 
Harper and workload ratings start out relatively flat and 
then degrade after a root location of 0.16 rad/sec. 
Comparison of Cooper-Harper and workload ratings for 


the height mode (fig. 10) shows similar results. The 
reduced Cooper-Harper rating of 7 at the height root loca- 
tion of 0.15 rad/sec was the only controllable run of the 7 
test runs and therefore does not represent a meaningful sta- 
tistical average. These results suggest that pilot workload 
was the primary consideration driving the Cooper-Harper 
rating. However, a question about the influence of pilot 
performance on Cooper-Harper ratings remains. 

Figure 11 shows the pilot’s performance, displayed as a 
percentage of desired performance, as a function of 
phugoid instability. The pilot performance on speed stays 
relatively flat until a root location of 0.3 rad/sec and then 
degrades, similar to the Cooper-Harper and workload rat- 
ings. However, the degradation in speed performance 
becomes no worse than approximately desired speed at the 
most unstable phugoid tested. The altitude performance 
data show some degradation as the instability increases but 
remains relatively flat and well within desired perfor- 
mance. Figure 9 shows even at the highest root locations 
of approximately 0.4 rad/sec, where control is marginal 
(Cooper-Harper ratings of 8-9), the performance data (fig. 
11) remain approximately less than or equal to desired per- 
formance. Therefore, despite some degradation in perfor- 
mance in speed, the pilots achieved levels of performance 
within the desired range even as the ratings deteriorated to 
the point of uncontrollability. 

Figure 12 shows the pilot’s performance, displayed as a 
percentage of desired performance, as a function of 
controllable height instability. Similar to the phugoid case, 
the pilots maintained desired performance on speed and 
altitude despite the degradation of Cooper-Harper ratings 
to 9 (fig. 10). 
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Figure 9. Average Cooper-Harper and workload ratings of the controllable cases 
as a function of phugoid instability (nontapping data). 
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Figure 10. Average Cooper-Harper and workload ratings 
of the controllable cases as a function of height instability (nontapping data). 
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Figure 11. Percentage of desired performance achieved by the pilot for nontapping, phugoid mode data. 



Figure 12. Percentage of desired performance achieved by the pilot for nontapping height mode data. 


The obvious conclusion is that pilot workload and not 
pilot performance is the major influence on Cooper- 
Harper rating. This result implies that a pilot flying aircraft 
with these unstable modes can achieve an acceptable level 
of performance, but the workload required to achieve this 
level of performance increases dramatically with instabil- 
ity. In measuring the flying qualities impact of these 
modes, a measure of workload is necessary. 

Figures 11 and 12 show that the speed performance was 
more critical than the altitude performance. With phugoid 


roots greater than 0.3 rad/sec, the speed performance 
degraded, whereas the altitude performance remained flat. 
Also, pilots achieved a lower percentage of altitude 
deviation than speed deviation across all instabilities. This 
lower percentage of altitude deviation may be influenced 
by the fact that equivalent airspeed is more difficult 
to control because it is a function of Mach and altitude. 
The choice of desired airspeed performance of ±5 knots 
equivalent airspeed (18 lbs/fit^ dynamic pressure) may also 
be an influence. Tight performance on speed was chosen 
because of the sensitivity of the propulsion systems to 
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dynamic pressure and because tight control of dynamic 
pressure is critical to the success of a singie-stage-to-orbit 
vehicle mission. 

Analysis of Pilot Technique for Large Instabilities 

In general, the data presented show that pilots can suc- 
cessfully control aircraft models with quite large phugoid 
and height mode instabilities. For example, phugoid insta- 
bilities with a time-to-double-amplitude of 1.7 sec were at 
times controllable, although the workload was very high. 
Inherent to the design of air-breathing hypersonic vehicles 
is the shaping of the aft fuselage to act as a nozzle. This 
design creates significant coupling between the vertical 
and longitudinal axes when the thrust is modulated. 
Figure 13 shows this principle with the linearized version 


of the GHAME at Mach 10. A positive amplitude throttle 
step produces a significant increase in vertical and longitu- 
dinal speed. This coupling characteristic is potentially 
quite useful to the pilot in controlling large instabilities. 
The following evaluation of the data determines whether 
or not the pilot took advantage of this coupling 
characteristic. 

Figure 14 shows the time history of a run split into three 
regions, A, B, and C, to show how one pilot controlled 
a phugoid instability with a time-to-double-amplitude of 
1.7 sec. At the beginning of region A, the pilot observed 
the negative trend on vertical speed and reacted with full 
positive aft stick deflection, thereby reducing the slope of 
the vertical speed but not negating it. At this point, the 
pilot would have lost control of the vehicle unless 
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Figure 13. Throttle step showing vertical speed and longitudinal speed coupling. 
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Region 



Time, sec 


Figure 14. Pilot technique controlling phugoid mode instability. 


something else was done to arrest the vertical speed. 
Region B shows the pilot’s reaction. The pilot increased 
the throttle sharply at the beginning of region B, thereby 
creating enough positive vertical speed to negate the 
downward trend. Once the pilot observed the sign reversal 
of the vertical speed, the pilot sharply decreased the throt- 
tle. 

At the end of region B, the pilot was confident control 
was retained and thus reduced the stick deflection off its 
limit. However, the vertical speed kept climbing positively 
and went unnoticed until approximately 160 sec, which is 
the beginning of region C. At this point, the pilot reversed 
the above procedure by reducing the stick deflection. 
Again, it was not enough to arrest the trend of the vertical 


speed, so the pilot applied a negative, sharp throttle incre- 
ment. Unfortunately, the throttle input was too late and the 
pilot lost control. 

Nevertheless, Figure 14 shows the pilot used 
sharp throttle inputs to take advantage of the coupling 
between longitudinal and vertical axes to control a 
vehicle otherwise uncontrollable. This technique with 
the throttle, a “bang-bang” technique, increased controlla- 
bility and minimized airspeed deviations. Note that 
the engine model for this simulator experiment was a 
simple, linear model with no time delay. A more 
realistic engine modeled with filters or time delays 
would reduce the effectiveness of this “bang-bang” 
technique. 


14 



Effect of Control Power on Flying Qualities 

Because the maximum lift coefficient was not exceeded, 
the inability to control vertical speed with full aft stick 
deflection (fig. 14) implies the pilot ran out of control 
power. However, the pilot supplemented the control power 
with the coupling of thrust between the vertical and longi- 
tudinal axes. The success of this technique implies the 
pilot was not limited by the capacity to handle the large 
instabilities, but rather by the available control power. If 
this implication is true, increasing the control power 
would allow the pilot to fly aircraft models with even 
higher instabilities. 

Figure 15 shows the impact of modifying the control 
power on one pilot’s Cooper-Harper ratings of phugoid 
instabilities. By varying the control surface deflection lim- 
its, three different control power variations are compared: 


the baseline control power, 0.5 of the baseline control 
power, and 1.5 of the baseline control power. The pilot 
controlled much higher levels of phugoid instability with 
the increased control power, as evidenced by the increase 
in the Level 2/3 border from root locations of 0.23 to 
0.32 rad/sec and Level 3/Uncontrollable border from root 
locations of 0.44 to 0.60 rad/sec. 

These data indicate that the long-term modes are flyable, 
even for large instabilities, as long as enough control 
power is available. As the available control power lessens, 
the bandwidth required increases and the workload 
necessary to maintain acceptable or controllable 
performance increases, resulting in an increase in 
Cooper-Harper and workload ratings. 

Figure 16 shows this relationship by plotting the vertical 
speed response of two phugoid instabilities to a step input 
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Figure 15. Effects of available control power on Cooper-Harper pilot ratings of phugoid mode instabilities. 
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Figure 16. Effects of control power on pilot workload requirements. 
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of stick deflection and the vertical speed control power limit, 
which is the point where full reversed stick deflection just 
arrests the vertical speed. Decreasing the time-to-double- 
amplitude reduces the reaction time necessary for the pilot to 
maintain controllability, thereby increasing the concentra- 
tion and effort required of the pilot. Decreasing the control 
power available, which is equivalent to lowering the “con- 
trol power limit” boundary (fig. 16), has the same effect. 
Therefore, the relationship between the unstable phugoid 
mode and the control power available impacts the flying 
qualities of a hypersonic vehicle and may influence the con- 
trol power requirements of the design. 

Not only could this hypothesis explain why the Cooper- 
Harper ratings were primarily a function of the workload 


rating, as discussed earlier; it also could explain the “cliff- 
like” ratings. For low levels of instabilities, the bandwidth 
required is low enough not to influence the pilot ratings. As 
the instability increases, or as the control power decreases, 
the bandwidth required begins to approach the maximum 
available from the pilot-vehicle system, resulting in deterio- 
rated ratings. 

Tapping Method Evaluation 

As stated earlier, the main objective of the tapping task 
was to measure the degree of PML, or workload, on the pilot 
as the instabilities increased. Since tapping was designed as 
a secondary task, no degradation in pilot ratings should 
occur in the main task because of tapping. Figures 17 and 18 
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Figure 17. Average Cooper-Harper ratings of tapping and nontapping data as a function of phugoid instability 
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Figure 18. Average Cooper-Harper ratings of tapping and nontapping data as a function of height mode instability. 
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compare tapping and nontapping Cooper-Harper ratings as 
a function of the phugoid and height mode instability lev- 
els respectively. Except for the rating at a real phugoid 
root of approximately 0.3 rad/sec (fig. 17), the ratings 
show little difference between tapping and nontapping 
data. From this observation, it can be concluded the sec- 
ondary task did not interfere with the main task. 

But the question on how well the tapping estimated pilot 
workload remains. Figure 19 shows the average of the 


pilot’s estimation of workload with tapping for 
the phugoid and height mode instabilities. As the phugoid 
instability increases, the workload stays relatively 
constant at first and then degrades sharply at a root loca- 
tion of approximately 0.3 rad/sec. The workload 
for increasing height mode instabilities exhibits a similar 
trend, degrading at a root location of 0.024 rad/sec. 
These observations match intuition as well, as one 
would expect workload to increase with increasing 
instability. 




(b) Height mode data. 

Figure 19. Average workload scale ratings as a function of height and phugoid instabilities (tapping data). 
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Figure 20 shows all PML data as a function of phugoid 
and height mode instabilities and overplots the average 
PML at each instability. The phugoid data (fig. 20(a)) 
shows an increase in PML after the root location of 0.3 
rad/sec, which was the point where the average pilot 
workload rating began degrading sharply. Nevertheless, 
the scatter in the PML data, which is as large as the 
previously noted increase in PML, obscures any trends 
exhibited by averaged PML data and prevents the 


conclusion that any significant correlation exists between 
increasing instability and PML. 

In addition, no correlation between PML and increasing 
instability is observed in the height mode data (fig. 20(b)) 
because the large degree of scatter in PML obscures any 
trends exhibited by the averaged PML data. Although 
consistent trends between PML and instability exist in 
some individual series of runs, these trends exhibit both 




(b) Height mode data. 

Figure 20. Perceptual motor load as a function of instability for height and phugoid mode data. 
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increasing and decreasing PML as a function of increasing 
instability. The lack of correlation between phugoid and 
height mode instabilities and PML is inconsistent with the 
pilot’s estimation of workload (fig. 19). The pilot’s estima- 
tion shows a strong correlation between workload and 
increasing instability. 

In some series of runs* PML increased as the number 
of the runs evaluated in the series increased. After the 
pilots evaluated a few runs in a series, they occasionally 
complained that their fingers became tired because of the 
tapping. This effect could increase the irregularity of their 
tapping and produce increasing PML with the number of 
the runs evaluated in the series. This factor may indicate 
pilot fatigue is sometimes a more significant factor in the 
tapping workload measurement than the increasing insta- 
bility. However, data from previous research 15 show no 
indication of degradation caused by fatigue in test runs of 
comparable length in time. No conclusive correlations 
were observed when the data were evaluated for individ- 
ual pilots. 

For this experiment, the tapping method was inconsis- 
tent in measuring pilot workload. No reason for the incon- 
sistency has been isolated to date because the various 
effects could not be separated out with certainty. 

SUMMARY 

The NASA Dryden Flight Research Center and the 
Technical University of Munich are cooperating to 
research the flying qualities impact of unstable long-term 
dynamics on planes in hypersonic flight. These issues 
were investigated with a dedicated flight simulator for 
hypersonic vehicles located at NASA Dryden. The 
results presented in this paper are as follows: 

1. Level 2 handling qualities could be maintained to a 
time-to-double-amplitude of 15 sec for the height 
mode and 3 sec for the phugoid mode. Level 2 is ap- 
propriate for the design of a manual backup control 
mode; however, the trends in the ratings indicated 
“cliff-like” tendencies that may warrant extra margin 
in design. The borders of controllability were a time- 
to-double-amplitude of 5 sec for the height mode and 
1.4 sec for the phugoid mode. A technique developed 
by the pilots to take advantage of the coupling of thrust 
between the vertical and longitudinal axes allowed 
these large instabilities to be flyable. 

2. The pilots were able to maintain a steady level turn at 
or within desired performance margins even at the 
highest levels of instabilities, but at a high cost of 
workload, as reflected by the increased Cooper-Harper 
ratings. 


3. The flying qualities impact from the long-term dynam- 
ics influence control power design requirements. The 
influence is shown by the deterioration of 
the Cooper-Harper ratings that accompanies decreas- 
ing available control power. As instability increases, 
the bandwidth required approaches the bandwidth the 
pilot is capable of producing with the aircraft. At this 
point, workload and Cooper-Harper ratings deteriorat- 
ed. Decreasing control power has a similar effect. 

4. A method for measuring pilot workload, the tapping 
method, did not work as shown by the inconsistent re- 
sults for these experiments. 
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Appendix 


LINEAR EQUATIONS AT MACH 10 AND AN 
ALTITUDE OF 110,000 FEET 

This appendix describes the linear equations of motion 
that were used in the flying qualities experiment. Basic 
aircraft equations and control system block diagrams are 
presented for both longitudinal and lateral-directional 
axes. 


LONGITUDINAL EQUATIONS 

Figure A-l shows the longitudinal block diagram. Pilot 
inputs into the control system are represented as longitudi- 
nal stick deflection, & gp , and throttle position, 8 r . The 
linear aircraft equations of motion are represented by the 
state-space equation, which is defined as follows: 
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Figure A-l. Longitudinal block diagram. 
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Fhugoid instabilities were introduced by varying the 
control system gain, K g1 with the values listed in table 2 
in order to obtain their corresponding instabilities. The 
rest of the control system gains remained constant. Height 
mode instabilities were introduced by varying the value of 
C D v (table 3). Varying C Q had the effect of varying only 
the element X in matrix A; all the other elements in both 
the A and B matrices were constant. Table A-l shows the 

^5 * ’ anc * ^ values for both phugoid and height mode 

variations. 

LATERAL-DIRECTIONAL EQUATIONS 

Figure A-2 shows the lateral-directional block diagram. 
Pilot inputs are described by 8 and 8 r for the lateral 
stick deflection and rudder pedal inputs respectively. The 
linear aircraft equations of motion are represented by the 
state-space equation. The lateral-directional control sys- 
tem gains and state-space matrices remained constant 
throughout both phugoid and height mode evaluations. 
The state-space matrices are as follows: 


Table A-l. Constants required for implementing height 
and phugoid mode variations. 



*8 

C D V 

X 


0 




4.2 



Phugoid mode variations 

10 

0 

-0.22xl0~ 2 


20 




25 




30 





0.00001 

- 0 . 89X10- 2 



- 0.00003 

0 . 18X10- 1 



- 0.00006 

0 . 38xl0 -> 

Height mode variations 

0 

- 0.00008 

0 . 51x10 -' 



- 0.00018 

0.12 



- 0.00024 

0.16 



- 0.00030 

0.2 



-0.79x10-’ 

0.19x10"’ 

0.96 

0.18X10" 5 

-0.42x1 Or 1 


-0.43x1 0“ 3 

-0.62x10-2 

3.76 

0.70X10" 5 

-0.16X10- 6 

A = 

0.143 

-0.99 

-0.17x10-' 

0.31x10-2 

0.30x1 0” 8 


1 

0.14 

oo 

X 

9 

-0.86x10-5 

0.61X10" 4 


0 

1 

-0.69x1 0" 4 

-0.53X10" 3 

0.86x10-5 


5.1 

1.2 


P 



0.20x10-’ 

-0.75x10"’ 


r 


B = 

0.91X10- 4 

0.86x1 0- 3 

X = 

p 

u = 


0 

0 


4> 



0 

0 


V 



57.3 

0 

0 

0 

0 


0 

57.3 

0 

0 

0 

C = 

0 

0 

57.3 

0 

0 


0 

0 

0 

57.3 

0 


0 

0 

0 

0 

57.3 


22 




950030 


Figure A-2. Lateral-directional block diagram. 
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